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ABSTRACT 

We have mapped the central region of NGC 4945 in the J = 2 — > 1 transition 
of 12 CO, 13 CO, and C 18 0, as well as the continuum at 1.3 mm, at an angular 
resolution of 5" x 3" with the Submillimeter Array. The relative proximity of 
NGC 4945 (distance of only 3.8 Mpc) permits a detailed study of the circumnu- 
clear molecular gas and dust in a galaxy exhibiting both an AGN (classified as 
a Seyfert 2) and a circumnuclear starburst in an inclined ring with radius ~2'.'5 
(~50 pc). We infer the systemic velocity ~585 km s" 1 from channel maps and 
PV-diagrams. We find that all three molecular lines trace an inclined rotating 
disk with major axis aligned with that of the starburst ring and large-scale galac- 
tic disk, and which exhibits solid-body rotation within a radius of ~5" (~95 pc). 
The rotation curve flattens beyond this radius, and the isovelocity contours ex- 
hibit an S-shaped asymmetry suggestive of a highly inclined bar as has been 
invoked to produce a similar asymmetry observed on larger scales. We infer an 
inclination for the nuclear disk of 62° ±2°, somewhat smaller than the inclination 
of the large-scale galactic disk of ~78°. The continuum emission at 1.3 mm also 
extends beyond the starburst ring, and is dominated by thermal emission from 
dust. If it traces the same dust emitting in the far-infrared, then the bulk of 
this dust must be heated by star-formation activity rather than the AGN. We 
discover a kinematically-decoupled component at the center of the disk with a 
radius smaller than 1'.'4 (27 pc), but which spans approximately the same range 
of velocities as the surrounding disk. This component has a higher density than 
its surroundings, and is a promising candidate for the circumnuclear molecular 
torus invoked by AGN unification models. 
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Subject headings: galaxies: active — galaxies: Seyfert — galaxies: starburst - 
galaxies: individual (NGC 4945) — galaxies: intergalactic medium — radio lines: 
galaxies 



INTRODUCTION 



NGC 4945 is a nearby, almost edge-on (i = 78°), disk galaxy (type SB(s)c d or SAB(s)cd) 
that is a member of the Centaurus Group of galaxies (IWebster et al.lll979l ). Its distance, 
determined most a ccurately from the lumino sity of stars at the tip of the red giant branch, 
is 3.82±0.31 Mpc (IKarachentsev et al.l 120071 ). This value is consistent with previous more 
cruder estimates of the distance to NGC 4945 (see the discussion in Bergman et al. 1992), 
and so henceforth we shall assume that 1" corresponds to 19 pc in NGC 4945. 

The central region of NGC 4945 contains an active galactic nuc leus (AGN) reyealed 



most unambiguously by its strong and variable hard X-ray emission (jlwasawa et al.l Il993 



Done et al.lll996t iMadejski et al.ll2000l ). H2O megamaser emission has been detected from 
the nucleus, likely distributed in a disk with radius of ~15 mas (~ 0.3 pc) around the 
AGN (Greenhill et al. 1997). Th e AGN is classified as a Seyfert 2 faraatz et all Il997i : 
Madejski et al.ll2000l : ISchurchll2002l ). consistent with fi nding s that Seyfert 2 but not Seyfert 1 
galaxies exhibit water megamasers (e.g., review by lLoll2005l ). Surrounding the AGN, there is 
an in clined circumnuclear starburst ring with a radius of ~2'/5 (~50 pc) seen most clearly in 
Pact (IMarconi et al.ll2000l ). The central region of NGC 4 945 is among the strongest and most 
prolific extragalactic sou r ces of molecular lines (e.g., iHenkel. Whiteoak. fc Mauersberger 
1994 ; ICurran et all l200ll ; IWang et all 120041 ) . For these reasons, NGC 4945 is a particu- 
larly attractive candidate for studying the nature of molecular gas at the center of an active 
galaxy, and the role this gas plays in fueling and perhaps also determining in part the ob- 
served properties of both the circumnuclear starburst and AGN. 

Because of its southerly declination (5 = —49°), studies of the molecular gas in NGC 4945 
have been largely restricted to si ngle-dish telescopes, in particular the 15-m Swedish-ESO 
Sub millimeter Telescope (SEST). IWhiteoak et al.l (jl990l ) made the first map of the central 
region of this galaxy in the J = 1 —> transition of 12 CO at angular resolution of 43" (mea- 
sured at full-width to half-maximum, FWHM, henceforth used throughout when quoting 
sizes). They attributed the observed strong central concentration of molecul ar gas to a disk 
with a radius of ~18" (340 pc). With just a single pointing at the nucleus, iBergman et al. 
( 119921 ) obtained high-quality spectra of both 12 CO and 13 CO in the J = 1 — > transition 



They found that the observed line ratio as a function of velocity can be modeled as a ho- 
mogenous ring in rigid-body rotation with an inner radius of ~8" (~150 pc) and an outer 
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radius of -15" (-280 pc). 



Dahlem et al.l (119931 ) mapped the central region of NGC 4945 in both the J = 1 — > 



and J = 2 — > 1 transitions of 12 CO at an angular resolution of 43" and 24" respectively. 
Based on the measured position- velocity diagram along the major axis of NGC 4945 in both 
transitions, they inferred the presence of a ring with a radius of 8"±3" (~150 pc) in the 
J = 2 —> 1 and 15"±3" (~280 pc) in the J = 1 —* transitions. They noted that this 
did not necessarily suggest the presence of two rings, but rather a radial variation in the 
excitation conditions and/or optical depth in a single ring with a given radial thickness. 



Mauersberger et al.l (119961 ) mapped the central region of NGC 4945 in the J = 3 



2 transition of 12 CO at an angular resolution of 15". In this transition, they found an 
intrinsic size for the central concentration of radius 5"2 ± 1"5 (~100 pc) along its major axis. 
Although the position- velocity diagram is once a gain consistent with a ring but with an even 
smaller radius than in lower 12 CO transitions, iMauersberger et al.l (119961 ) were careful to 
point out that other structures (e.g., disk, bar, spiral, or even two independent molecular 
concentrations) are possible. 

All the abovementioned single-dish observations infer the presence of a molecular ring 
with size smaller than the respective angular resolutions of these observations. It is not clear 
how this molecular ring is related to the starburst ring, whether it comprises the reservoir for 
fueling the active nucleus, and whether it is relate d to the hypothe tical circumnuclear molec- 
ular torus required by AGN unification models (lAntonuccil 1 19931 ) . Observations at higher 
angular resolutions are needed to spatially resolve and hence clarify the spatial-kinematic 
distribution of the mole cular gas at the center of NGC 49 45. The first such observations were 
made quite recently by 1 Cunningham fc Whiteoakl (120051 ). who used the Australia Telescope 
Compact Array (ATCA) to map a central region of radius 17" (323 pc) in the J = 1 —> 
transitions of HCN, HCO + , and HNC at an angular resolution of 5"6 x 3"5 (106 x 66 pc). 
These lines trace molecular gas at about an order of magnitude higher density (~10 4-5 cm -3 ) 
than the abovementioned 12 CO lines (~10 3 ~ 4 cm -3 ). The HNC map, which best traces the 
overall spatial-kinematic distribution of the molecular gas, reveals an inclined rotating disk- 
like feature with a radius of 4'.'25 ± 0"25 (~80 pc) along the major axis and position angle 
of 64°, somewhat misaligned from the major axis of both the starburst ring and large-scale 
galactic disk at a position angle of ~45°. Based on the measured position- velocity diagram, 
Cunningham fc Whiteoakl (120051 ) argue that the HNC line traces a ring with an inner radius 
of 3" (57 pc). Their maps, however, are affected by a lack of short baselines that result 
in an inability to detect structures larger than ~15" (~285 pc). Furthermore, line absorp- 
tion against the strong central continuum source at their observing wavelength of 3.3 mm 
(90 GHz) may compromise the detectability of the innermost features. 
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Here, we present the first interferometric observations of the central region of NGC 4945 
in 12 CO, as well as in 13 C0 and C 18 0, at their J = 2 — > 1 transitions. The angular resolution 
attained is much higher than previous single-dish CO observations, and we properly resolve 
for the first time the central molecular gas concentration in CO lines. Our maps recover 
virtually all the emission detected in single-dish observations, and are not affected by line 
absorption against the central continuum source. We show that the CO gas is distributed 
in an inclined disk exhibiting rigid-body rotation out to a radius of ~5" (~95 pc); i.e., the 
molecular disk extends beyond the circumnuclear starburst ring seen in Paa. The rotation 
curve flattens beyond this radius, and the isovelocity contours exhibit an S-shaped asymme- 
try suggestive of a highly inclined bar as has been invoked to produce a similar asymmetry in 
both molecular and atomic hydrogen gas extending to much larger radii. We do not detect 
any central hole in the disk, placing an upper limit of 2".2 (42 pc) on the radius of any such 
hole. Instead, we detect for the first time a spatially-unresolved component at the center of 
the disk that is kinematically decoupled from the surrounding disk. This inner component 
exhibits a broad velocity width comparable with the overall range in velocities exhibited by 
the surrounding disk, and is a good candid ate for the hypot hesized circumnuclear molecular 
torus invoked by AGN unification models (j Antonuccil 1 1 9931 ) . 



2. OBSERVATIONS AND DATA REDUCTION 



We observed NGC 4945 with the Submillimeter Array (SMA; [Ho et all 120041 ) . located 
on Mauna Kea in Hawaii, on 2004 May 22 and 24. To achieve the most circular beam 
possible, we used the compact configuration of the array that is elongated in the north-south 
direction tailored for southern hemisphere sources. Over the course of the observations, this 
array configuration provided projected baselines ranging from about 6 m to 74 m. Seven 
of the eight antennas of the array were in operation at the time of our observations. Each 
of the SMA antennas has a diameter of 6 m, and at our observing wavelength of 1.3 mm 
(frequency of 230 GHz) a corresponding FWHM primary beam of 52" (i.e., field of view of 
nearly 1 kpc). 

The array was pointed toward the center of NGC 4945 as defined by the position of H2O 
megamasers at right a scension a(2000) = 13 h 05 m 27:5 and declination 5(2000) = 49°28' 05"56 



(IGreenhill et al.lll997f ). We observed the 12 CO, 13 CO, and C 18 lines in their J = 2 -> 1 
transitions simultaneously by placing the 12 CO line in the upper sideband (USB) and the 
13 CO and C 18 lines in the lower sideband (LSB) of the double-sideband receivers. The 
correlator was configured to provide a total bandwidth of 2 GHz in each sideband using 
twenty-four independent spectral windows ("chunks"). Each window is divided into 128 
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channels with a channel width of 0.8125 MHz (i.e., ~1.1 km/s), giving a total frequency 
coverage in each sideband (with a small overlap between windows at their edges) of 1.9 GHz 
(i.e., -2570 km/s). 

Because of its low declination, NGC 4945 is visible for only ~2 hrs each day from the 
SMA, and never exceeds an elevation of ~25°. Although the weather was good (zenith 
opacity r = 0.15 — 0.20) on both days of our observations, because of its low elevation 
the sky opacity towards the source was relatively high resulting in system temperatures of 
400 K to 600 K (single sideband). On each day, we observed NGC 4945 for a total on- 
source integration time of 53 mins. We used two quasars for complex gain (i.e., amplitude 
and phase) calibration. The stronger quasar 1334 — 127, which lies at a larger angular 
separation of ~40° from NGC 4945, was used for amplitude calibration. The weaker quasar 
1313 — 333, which lies at a smaller angular separation of ~20° from NGC 4945, was used for 
phase calibration. The flux density of the amplitude calibrator 1334 — 127 was scaled from 
observations of Callisto. The quasar 3C279 was observed for bandpass calibration. 

We calibrated the data in the standard fashion using the software package MIR-IDL, 
and made images using AIPS. The results shown in this paper correspond to the data taken 
on May 24 only. Data taken on May 22, at a different sidereal time, have a different uv- 
coverage (greater number of longer projected baselines) and deep negative sidelobes that 
made imaging of the line emission difficult. To make a map of the continuum emission, we 
used the line-free channels (total bandwidth of 3 GHz) in both the LSB and USB. To make 
maps in line emission, we first used the line-free channels in the visibility plane to derive 
the continuum emission, which was then subtracted from all the channels. The synthesized 
beam achieved with natural weighting of the visibilities is 5'/l x 2'/8 (~100 x 50 pc) at a 
position angle of 8?8. 



3. RESULTS AND INTERPRETATION 



3.1. Continuum Emission 



3.1.1. Contribution from dust 



In Figured], we show a map of the 1.3-mm con tinuum emission supe rposed on a Pact 
image of the central region of NGC 4945 taken from lMarconi et al.l (120001 ) . The continuum 
source has a total flux density of 1.3 ± 0.2 Jy, and is clearly resolved. Gaussian fitting 
yields a deconvolved size (at FWHM) of 9?8 x 5'.'0 (±0'/3) (186 x 95 pc) with the major axis 
at a position angle (PA) of 28° ± 3°. The centroid of this source is located at a(2000) = 
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13 h 05 m 27!59 ± 0!01 and 5(2000) = -49°2 8' 06? 1 ± 0"2, which is l'/8 to the south-east of the 
nominal centroid of the H 2 megamasers (jGreenhill et al.ll 19971 ). The latter coincides within 
measurement uncertainties with the centroid of the central radio continuum source detected 



at cm- wavelengths (jElmouttie et al.lll997l ). and presumably marks the location of the AGN. 



The major axis of the 1.3-mm continuum source is not aligned with the major axis of 
either the starburst ring or larger-scale galactic disk, both of which have a PA « 45°. Along 
the major axis of the Paa starburst ring, the 1.3-mm continuum extends beyond the inner 
bright rim of the starburst ring at a radius of ~2'/5 (~50 pc), and also beyond the detectable 
outer radius of this ring at ~5" (~100 pc). On the south-eastern side of the ring, the 1.3-mm 
continuum clearly extends beyond the measured extent of the Paa emission. 

The central region of NGC 4945 has been imaged in the continuum at centimeter wave- 
lengths at angular resolutions comparable with that attained here. At 21 cm (1.4 GHz), the 
source has a deconvolved siz e of 7 / . / 6 x 3'.'4 (± 0?2) with major axis at PA = 42°. 5, and a total 
flux density of 4.6 ±0.1 Jy (jOtt et al.ll200ll ). At 5 cm (6 GHz), the source has a smaller de- 
convolved size of 5" 7 x 2!0 (±ff!l) at PA = 43° ± 1°, and a total flux density of 2.04 ±0.04 Jy 
(IWhiteoak fc Wilsonlll990i ). The size of the continuum source at centimeter wavelengths is 
therefore much smaller than that measured at 1.3 mm. Instead, the centimeter continuum 
source has a size comparable with the starburst ring, and its major axis is aligned with that 
of the starburst ring. The steep negative spectral index of this source indicates that non- 
thermal (synchrotron) emission dominates at centimeter wavelengths, and presumably arises 
from star- formation- related activity (e.g., supernovae) in the starburst ring (in addition to 
any unresolved emission from the AGN). 



Cunningham fc Whiteoakl (120051 ) have imaged the central continuum source at 3.3 mm 



(90 GHz) at an angular resolution comparable with that attained here. Like us, they find 
an elongated source whose centroid is offset to the south-east of the H 2 megamasers, and 
whose major axis is at PA « 29°. With a reported deconvolved source size of 7'.' 6 x 2'.'0 and 
peak flux density of 0.13 Jy, the corresponding total flux density for a Gaussian source is 
~1.0 Jy. Extrapolating from 21 cm and 5 cm, the expected flux density at 3.5 mm is ~0.5 Jy, 
only half that actually measured. Extrapolating to 1.3 mm, the estimated contribution from 
nonthermal emission is ~0.3 Jy, less than one quarter the flux density measured at this 
wavelength. Even if we assume that nonthermal emission dominates at 3.3 mm, this emission 
would only contribute at most about half the total flux density measured at 1.3 mm. Thus, 
we conclude that dust emission dominates at 1.3 mm (and also contributes significantly at 
3.3 mm), and that this explains the different dimensions and position angles of the source at 
millimeter and centimeter wavelengths. In §4.1[ we discuss the implications of these results 
for the origin of the central far-IR emission from NGC 4945. 
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3.1.2. Molecular gas mass from dust 

To estimate the mass of molecular gas from the inferred dust emission at 1.3 mm (assum- 
ing a gas to dust ratio of 100), we first subtract the estimated cotribution from nonthermal 
emission (~0.3 Jy) from the total continuum emission (~1.3 Jy). The gas mass is then given 
by 



MM = 17 30 ^gi»I^ 
T d [K\K d {p)[ g 1 cm 2 J 



(1) 



(]Hildebrandlll983l ) . where the dust continuum flux density S u ~ 1 Jy, distance D = 3.82 Mpc, 



wavelength A = 1.3 mm, and dust absorption coefficient Kd(1.3mm) w 3 x 10 3 g 1 cm 2 . As 



suming a dust temperature T d « 40 K as inferred from far-infrared measurements (IBrock et al 
19881 ). we find that M gas w 3.6 x 10 8 M Q . If we adopt a dust temperature corresponding 



to the peak brightness temperature measured in 12 CO(2-l) of ~30 K ( §3.21) . the gas mass is 
then M gas « 4.7 x 10 8 M . As we shall show in §3.2.31 the gas mass inferred from dust is 
comparable with values inferred from the three observed CO lines. 



3.2. Line emission 

3.2.1. Spatial-kinematic structure 

Like the continuum emission, the 12 CO(2-l), 13 CO(2-l), and C 18 0(2-l) lines are strongly 
concentrated towards the center of NGC 4945. Figures [2] and [3] show channel maps of the 
12 CO(2-l) and 13 CO(2-l) emission smoothed to a velocity resolution of ~20 km s _1 , with 
the circumnuclear starburst ring seen in Paa plotted as an ellipse. In Figure HI we show 
the corresponding total intensity (moment 0) as well as intensity-weighted mean-velocity 
(moment 1) maps. 

An inspection of these maps reveals that the molecular emission originates (primarily) 
from a highly-inclined rotating disk extending beyond the Pact starburst ring. Gaussian 
fitting to the moment maps yields a deconvolved size (at FWHM) for the disk of 16'.'4 x 10'.'8 
(±0"1) (~310 x 205 pc) and major axis at PA= 38° ± 1° in 12 CO(2-l), a size of 14"2 x 6"6 
(±0"1) (~270 x 125 pc) and PA= 43° ± 1° in 13 CO(2-l), and a size of 13"1 x All (±0'.'l) 
(~250 x 90 pc) and PA= 45° ± 1° in C 18 0(2-l). The derived position angle of the major 
axis in the more optically thin 13 CO(2-l) and C 18 0(2-l) lines is in good agreement with 
the position angle along which the velocity gradient of the disk reaches a maximum in these 
lines. The major axis of the disk (mean PA = 44°) is therefore well aligned with that of the 
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starburst ring as well as the larger-scale galactic disk at a position angle of ~45°. We see no 
evidence for the central hole in our data, suggesting that any central hole in this disk has 
a size smaller than the synthesized beam, which measures 4'.'4 along the major axis of the 
disk. 

In all three molecular lines, the inner region of the disk within a radius of ~5" (~95 pc) 
appears to exhibit simple circular rotation with the isovelocity contours perpendicular to 
major axis of the disk as can be seen in Figure HI Beyond this radius, however, the disk 
exhibits significant deviations from circular rotation as can be best seen in 12 CO(2-l). Here, 
the isovelocity contours on the north-eastern side of the disk twist to the north, and on 
the south-western side twist to the south, producing a S-shaped asymmetry In this outer 
region, the spatial structure of the disk closely resembles the distribution of the 1.3- mm 
continuum emission (e.g., the south-east extension), reinforcing our earlier argument that 
the continuum emission at 1.3 mm is likely dominated by dust. 

The position-velocity (PV-) diagrams of all three molecular lines along the inferred 
major axis of the disk are shown in Figure The PV-diagrams along the major axis 
in 13 CO(2-l) and C 18 0(2-l) indicate (primarily) rigid-body rotation, while in 12 CO(2-l) 
it shows additional complex emission with a broad velocity width at or near disk center. 
Cuts along different position angles reveal that this "excess" emission is always present as a 
(vertical) strip in velocity at the origin of the PV-diagram, not just in 12 CO(2-l) but also 
in 13 CO(2-l) and C 18 0(2-l). For example, in FigureOwe also show the PV-diagrams along 
the minor axis of the disk (PA = 134°), as well as halfway between the major and minor 
axes of the disk (PA = 89°), in all three molecular lines. Armed with this knowledge, a 
closer inspection reveals that this feature can just be seen in the PV-diagrams along the 
minor axis of the disk in 13 CO(2-l) and C 18 0(2-l). This feature therefore corresponds to 
a kinematically-decoupled component with a size smaller than the synthesized beam (i.e., 
projected radius as small as l'/4 or 27 pc) located at the center of the disk. In the channel 
maps, this kinematically-decoupled feature appears as the ever-present emission towards the 
center of the disk even at the most blueshifted and redshifted velocities detectable. By 
contrast, in an inclined disk exhibiting only rigid-body rotation, the most blueshifted and 
redshifted velocities should originate from just the outermost regions of the disk along its 
major axis. 

We infer from both the channel maps and PV-diagrams a systemic velocity of ~585 km s _1 
measured with respect to the local standard of rest. (All velocities quoted here are relative to 
the local standard of rest, which is 4.6 km s -1 lower than the heliocentric velocity.) By con- 
trast, single-dish observati ons in 12 CO(1-0) and 12 CO(2-l) infer a significantly lower systemi c 



velocity of ~558 km s" 1 flWhiteoak et al.l Il990l : Bergman et al.l Il992l ; iDahlem et al.lll993f ). 
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The HNC observation of ICunningham fc Whiteoakl (120051 ) indicate a systemic velocity, as 
measured at the midpoint between the two strongest HNC peaks in the PV-diagram (their 
Fig. 11), of ~570 km s _1 , about halfway between single-dish and our interferometric CO 
measurements. Th e systemic veloc ity inferred from interferometric observations in atomic 
hydrogen (HI) gas (jOtt et al.ll200ll ). averaged over the entire galaxy, is 557 ± 3 km s" 1 , in 
agreement with the abovementioned value derived from single-dish CO measurements. On 
the other hand, HI detected in absorption toward s the centroid o f the central continuum 
source has a velocity of ~585 km s _1 (see Fig. 5 of lOtt et al.ll200ll ). similar to the systemic 
yelocity inferred here. Arcsecond imaging of hydrogen recombination lines near 8.6 GHz 
(IRoy et al.ll2005l ) shows that the emission peaks at the location of the central radio contin- 
uum source; these recombination lines have a systemic velocity of ~576 km s _1 , again close 
to that inferred here. We therefore conclude, based on observations that better trace gas 
closer to the center of the galaxy, that the systemic velocity of NGC 4945 is more likely 
about 585 km s" 1 . 



3.2.2. Comparison with single-dish observations 



The CO(2-l) emission towards the center of NGC 4945 has been observed a num- 



ber of times with the SEST ( Dahlem et al.1 Il993l ; iHenkel. Whiteoak. fc Mauersbergerl 11994 



Curran et al.ll200ll ; IWang et al.1 12004 ) . These observations often give different line profiles 
and intensities, caused most likely by inacc u racies in te lescope pointing. Only the 12 CO(2- 
1) lin e profiles measured by iDahlem et al.1 (119931 ) and IHenkel. Whiteoak. fc Mauersberger 
( 119941 ) appear similar, and which also agree with that measured here. In Figure El we plot 
our spatially-integrated 12 CO(2-l) line profile corrected for the primary beam of the SMA 
and convolved to the primary beam of SE ST, together with the line profile measured by 
Henkel. Whiteoak. Mauersbergerl (jl_994j). The integrated line intensity (in main beam 



brightness temperature) we measure of 9 25 ± 13 K km s -1 is ~90% of that measured by 
Henkel. Whiteoak. fc Mauersbergerl (119941 ) of 1050 ±6 K km s _1 (as quoted in IWang et al. 
20041 ) . The quoted uncertainties do not take into account any errors in flux calibration, and 
so we have likely recovered the bulk if not all of the 12 CO(2-l) emission present in the same 
region. 

The 13 CO(2-l) and C 18 0(2-l) emissions to ward the center of NGC 4945 also have 



Curran et al. 


2001; 


Wang et al. 


2004) 



'he 13 CO(2-l) line profile we measure is similar in 
shape to that measured i n 12 CO(2-l) as shown in Figure El and also simil ar to the 13 CO(2-l) 



line profiles measured by IHenkel. Whiteoak. &: Mauersbergerl (119941 ) and IWang et al.l (120041 ) 



- 11 - 



but not lCurran et al.l (120011 ) . The integrated lin e intensity we measure of 82 ±3 K km s 1 is 



similar to that measured by Wang et al. (ho04) of 81.2 ± 0.7 K km s 1 (wh ich has a higher 
S/N ratio than that measured by lHenkel. Whiteoak. fc Mauersbergerlll994j ). 



The C 18 0(2-l) line profile we measure is again si milar in shape to th at measured in 
12 CO(2-l), but is different to those measured by either ICurran et al.l (120011) or IWang et al. 



(2004). The line profile we measure is closest in shape to that measured by lHenkel. Whiteoak. &: Mauersberg 



(jl994l ). The integrated line intensity we measure of 31 ± 3 K km s 1 also is similar to that 
measured by lHenkel. Whiteoak. &: Mauersbergerl (119941 ) of 32 ± 2 K km s -1 . In Table 1, we 



summarize the integ rated intensities that we measure fo r all t hr ee lines toget h er wit h the 
values measured b y iHenkel. Whiteoak. fc Mauersbergerl (119941 ). ICurran et al.l (120011 ). and 
Wang et all (12004 ) with the SEST. 



3.2.3. Molecular gas mass from CO lines 



To estimate the mass of molecular gas traced in 12 CO(2-l), we first use the conversion 
factor between the brightness temperature of the 12 CO(1-0) line and molecular hydrogen 
column density of 1 x 10 20 (K km s -1 )~ ^ cm" 2 as has been p roposed to be appropriate in 
the central regions of galaxies (e.g., see iPaglione et al.ll200ll ). Table lists the brightness 
temperature of the 12 CO(2-l) line if its emission as mapped with the SMA had been observed 
by SEST (details in §3.2.21) . In Table [5], we list the actual measured brightness of the three 
observed lines. Assuming a line ratio (in main beam brightness temperature) of 12 CO(2-l) 



to 12 CO(1-0) of 1.2 ± 0.1 (from lDahlem et al.lll993l . after correcting for the different beam 



sizes), we derive a molecular hydrogen column density of 5.9 ± 0.1 x 10 22 cm -2 . Given the 
measured source size as described in §3. 2. H and listed also in Table the corresponding mass 
in molecular hydrogen gas is therefore (1.63 ± 0.03) x 10 8 M Q . 

To compute the molecular gas mass from the 13 CO(2-l) and C 18 0(2-l) lines, we as- 
sume that this gas is in local thermal equilibrium (LTE). We adopt the abundanc e ra- 



tio [ 12 CQ1/[ 1 3 CQ1= 50 and [ 12 CO]/[C 18 0]= 200 as estimated bv ICurran et al.l fl200l|) and 



Wang et al.l (12004 ). We assume an excitation temperature of ~30 K based on the peak 
brightness temperature measured for the 12 CO(2-l) line (i.e., 16.5 Jy/Beam), and that the 
13 CO(2-l) and C 18 0(2-l) lines are optically thin. In this way, we derive a column density 
from the 13 CO(2-l) line of (6.7 ± 0.3) x 10 22 cm" 2 , and column density from the C 18 0(2-l) 
line of (1.65 ± 0.17) x 10 23 cm~ 2 . Their corresponding molecular hydrogen gas masses are 



(1.03 ±0.05) x 10 8 M as traced in 13 CO(2-l) and (1.21 ±0.1) x 10 8 M as traced in C 18 0(2- 
1). The estimated mass in molecular gas from all three lines agrees with each other, and 
are also roughly comparable with the value estimated from the dust continuum emission of 
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-4.7 x 10 8 M (§3X2]). 

We also compute the approximate mass in molecular hydrogen gas within the region 
where the disk exhibits circular rotation (i.e., radius < 5"). Table [5] lists the measured 
brightness temperature within a Gaussian region of size 10" x 4.7" and position angle of 
44 centered on the disk. The inferred column density and corresponding mass in molecular 
hydrogen gas are ~ 10 23 cm -2 and ~ 10 7 M as listed in Table [5j Follow the same method 
mentioned above, and the emitting area of 10"by 4"7 (assume the inclination of the disk is 
62°). The gas mass derived from 12 CO is (3.22 ± 0.06) x 10 7 M , and the gas mass derived 
from 13 CO and C 18 is (3.0 ± 0.1) x 10 7 M and (3.1 ± 0.5) x 10 7 M , respectively. 



3.2.4- Density and Temperature 



To search for any radial variations in density and temperature of the molecular gas, we 
computed line ratios (in measured brightness temperature) in 12 CO/ 13 CO, 12 CO/C 18 0, and 
13 CO/C 18 along the major axis of the disk. The results are shown in Figured Both the 
12 CO/ 13 CO and 12 CO/C 18 line ratios increase from the north-eastern (redshifted) side to 
the south-western (blueshifted) side of the disk, spanning the range 5-15 in 12 CO/ 13 CO and 
13-30 in 12 CO/C 18 0. On the other hand, the 13 CO/C 18 ratio remains roughly constant 
along the major axis of the disk, occupying a relatively narrow range between 2.2 and 3.0. 

To determine the line ratios of the inner spatially-unresolved but kinematically-decoupled 
component described in §3.2. H we turn to the PV-diagram along the minor axis of the disk as 
shown in Figure [SI Notice that the emission from this inner component away from the veloc- 
ity range 520-650 km s -1 , which straddles the systemic velocity, can be well separated from 
that of the surrounding disk. In this way, we compute for the inner component line ratios 
that span the range 12-17 in 12 CO/ 13 CO, 17-26 in 12 CO/C 18 0, and 1.2-1.5 in 13 CO/C 18 0. 
By comparison, the corresponding line ratios along the minor axis at velocities in the range 
520-650 km s -1 , where emission from the inner component and surrounding disk cannot be 
easily separated, are generally lower, spanning the range 8-15 for 12 CO/ 13 CO, 14-18 for 
12 CO/C 18 0, and 1.2-1.6 for 13 CO/C 18 0. 



We have used the LVG approximation (IGoldreich fc Kwanlll974l ) to compute the physical 
conditions of the molecular gas implied by the measured line ratios. The collision rates for 
CO in the te mperature ra n ge 10 -250 K were taken from lFlower fc Launayl (119851 ) and 500- 



2000 K from lMcKee et a 

[ 13 co]/[h 2 ; 



rl2 



CO]/[C 18 0] 



Jl982h. In these ca lculations, we adopted a relative abundance of 



1 x 10 6 (Solomon et a 



197a ). and isotopic ratios [ 12 CO]/[ 13 CO] = 50 and 



200 JWang et alil2004h . We assume that the emission in all lines are emitted 



13 



from the same region (i.e., a one- zone model), and a velocity gradient dv/dr w 1 km s _1 pc -1 
given the measured linewidth and diameter for the 13 CO emitting region of ~285 km s _1 
and ~270 pc respectively. The results for all three sets of line ratios are shown together in 
Figure [HI and exhibit the following trends: (1) at temperatures T < 100 K, all three line 
ratios exhibit only a weak dependence on temperature, with lower ratios corresponding to 
higher densities; (2) at temperatures T > 100 K, for a given line ratio the density increases 
with temperature, with smaller line ratios continuing to indicate higher densities; and (3) the 
line ratios 12 CO/C 18 and 13 CO/C 18 tend to indicate a higher density than 12 CO/ 13 CO. 

The measured line ratios therefore imply a decrease in gas density and column den- 
sity from the north-eastern (redshifted) side to the south-western (blueshifted) side of the 
disk. Our measurements do not place strong constraints on the gas temperature, but 
we note that the measured 12 CO/ 13 C O line ratios are comparable with average values of 
~13 ± 5 found in starburst gal axies (lAalto et al.l 119951 ) and ~13 ± 1 in Seyfert galaxies 
(IPapadopoulos fc Seaquistlll998l ). If the similar line ratios indicate similar physical condi- 
tions, then the preferred solution is for temperatures T < 100 K and densities n(fi2) ~ 
10 3 cm" 3 , typical values of the bulk properties for giant molecular clouds in our Galaxy. In 
this regime, the gas density changes by a factor of just a few across the major axis of the disk. 
Consistent with this idea, single-dish observations in mul tiple lines and transitions find a den- 
sity of ^(3-5) xlO 3 cm" 3 and temperature of ~100 K (IHenkel. Whiteoak. fc Mauersberger 



1994 ; ICurran et al.ll200ll ; IWang et al.ll2004l ) for the bulk of the centrally unresolved molecu- 
lar gas. At radii beyond ~10" (~190 pc) at or beyond the outer regions of the disk, where 
13 CO and C 18 are not detectable at our sensitivity limits, the line ratios 12 CO/ 13 CO and 
12 CO/C 18 have lower limits that are larger than their measured values at smaller radii. 
This indicates that the CO emitting outer region of the disk is dominated by relatively 
diffuse gas at densities of ~10 2 cm -3 . 

The inner kinematically-decoupled component exhibits a 13 CO/C 18 line ratio that 
is near unity (1.2-1.5); i.e., the less opaque C 18 line has nearly the same flux density 
as the 13 CO line. The 12 CO/C 18 line ratio of 14-18 also is relative ly low compared 



Sage et al 



1991 



with values in the range 40 -70 typically found in other galaxies (see 
Henkel &: Mauersbergerl Il993l ). At face value, this suggests that even the 13 CO and C 18 
lines are nearly optically thick in the inner kinematically-decoupled component. Alterna- 
tively, the C 18 species may be unusually abundant within the inner component. 

We consider the second possibility first. 18 enrichment can occur through its produc- 
tion from 14 N by He-burning in high-mass stars ( > 8 M n ), and subsequently dispersed in 
Wolf-Rayet phase or in ty pe II supernova explosions ( Sage et aDl991 ; Henkel fc Mauersberger 
19931 ; lAmari et al.l 119951 ). Such enrichment may have occurred during past circumnuclear 
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starsburt activity in NGC 4945, analogous to that currently seen as the P a ring. Such low 
line r atios are seen in at least one other star forming galaxy, NGC 6946 (IMeier fc Turner 



20041 ) . On the other hand, not all star for ming galaxies show su ch low line ratios; e.g., 
both M82 JWeiss et alJboOlh and NGC 253 JSakamoto et al.ll2006h exhibit 12 CO/C 18 or 
CO/C 18 line ratios comparable with the mean observed in other galaxies. 



13/ 



We turn back to the first possibility, which requires both the 13 CO and C 18 lines to be 
optically thick. To infer the required physical properties of the molecular gas, we use the same 
LVG approximation and the same assumptions as before for the surrounding disk, except 
that we now assume an order of magnitude higher velocity gradient as is more appropriate 
for the inner component. The results are shown in Figure [TUl which reveals that the physical 
properties of the gas derived from the 12 CO/ 13 CO and 12 CO/C 18 line ratios are comparable 
but very different from those derived from the 13 CO/C 18 line ratio. Our results imply that 
the one-zone model is not valid for the inner kinematically-decoupled component, and that 
the 13 CO(2-l) and C 18 0(2-l) lines trace a different denser region than the 12 CO(2-l) line. 
This situation is not uncommon in galaxies, with the 12 CO emission originating from more 
exten ded and diffuse gas and the 13 CO or C 18 Q emission fro m more compact and dense gas 



[e.g.. iDownes et al.lll992t IWall et al.lll993l ; lAalto et al.lll995l ). The physical properties of the 
molecular gas as inferred from the 12 CO/ 13 CO and 12 CO/C 18 line ratios are comparable 
with those inferred in the surrounding disk, and presumably correspond to the more diffuse 
part of the inner component. On the other hand, the measured 13 CO/C 18 line ratio implies 
that even at temperatures as low as T « 10 K, the gas density is ~5 x 10 4 — 1 x 10 5 cm -3 . 
This is 1-2 orders of magnitude higher than the gas density in the surrounding disk, and 
presumably corresponds to the denser part of the inner component. The inferred gas density 
of the denser part increases towards higher temperatures, with densities about an order of 
magnitude higher still at T « 100 K. 



4. DISCUSSION 



4.1. Dust Heating 



NGC 4945 is one of the three brightest IRA S point sources b eyond the Magellanic 
clouds. It has a far-IR luminosity of ~2 x 10 10 L Q ( jBrock et al.lll988l ). which is comparable 
with that radiated at all other wavelengths combined. Nearly all (at least 80%) of the far-IR 
emission arises from a central region no larger than 12" x 9" (230 pc x 170 pc), which is 
comparable in size to the central continuum source that we detected at 1.3 mm. The far-IR 
emission from this central source is attributed to dust at a temperature of ~40 K. 
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As pointed out by lMarconi et al.l (120001 ). it is not clear whether the central dust-emitting 
region in NGC 4945 is heated by the AGN or circumnuclear starburst. This situation reflects 
the general difficulty in deducing the nature of the source that heats dust in the nuclear region 
of active galaxies. In the case of NGC 4945, our observation reveals that the central dust that 
emits at 1.3 mm (size ~190 pc along the major axis) spans the entire observable extent and 
somewhat beyond the circumnuclear starburst ring. This spatially-extended dust is unlikely 
to be heated predominantly by the AGN as the latter is embedded in obscuring material 
that prevents UV photons (13.6 < hy < 500 ev) from p enetrating beyond a distance of at 
most 1'.'5 (~30 pc) from the center (jMarconi et al.ll2000l ). Even soft X-rays from the AGN 
can only escape along a X-r ay plume (beli eved to be blown by a nuclear starburst) that 
emerges north-west of center (jSchurchll2002l ); i.e., the axis of the X-ray plume is orthogonal 
to the plane of the central molecular gas and dust disk . Thus, if the dust emitting at 
1.3 mm can be used as a proxy for that emitting at far-IR wavelengths, the bulk of this dust 
is likely heated by the circumnuclear starburst. Observations at higher angular resolutions 
are required to search for any dust heated by the central AGN and to study the very inner 
structure. 



4.2. Central Molecular Concentration 

4-2.1. The central disk 

Our observation spatially resolves the central molecular gas concentration as traced in 
CO into an inclined rotating disk. The radius of this disk, measured over the region where 
it exhibits rigid-body rotation, is ~5" (~95 pc), although the central CO-emitting region 
extends beyond this radius. The overall radial size of the emitting region in 12 CO(2-l) is 8'.'2 
(156 pc), which i s simil ar to that inferred from spatially-unresolved single-dish observations 



by lDahlem et al.l (119931 ) . Their PV-diagram along the major axis of NGC 4945 (their Fig. 3) 
shows two local intensity peaks at velocities of ~430 km s -1 and ~710 km s _1 separated by 
16" ± 6", interpreted as the two cross-sections of a highly-inclined ring. Our PV-diagram 
(Fig. [5]) also shows two local intensity peaks near these velocities, specifically at 470 km s _1 
and 730 km s _1 , separated (as measured from their centroids) by ~13". These local intensity 
peaks in the PV-diagram correspond to the location where the rotation curve changes from 
rigid body to nearly flat, rather than tracing local spatial peaks corresponding to the two 
cross-sections of an inclined ring. Our spatially-resolved observations show no evidence for a 
central hole (or depression) in 12 CO(2-l), nor in 13 CO(2-l) and C 18 0(2-l), thus placing an 
upper limit of 2'/2 (42 pc) on the radius of any such hole. 

The central molecular concentration is therefore a disk or torus (if it has a spatially- 
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unresolved central hole) rather than a ring. If it has a thickness much smaller than its 
observed dimensions, then the disk must have an inclination of 62° ±2° to the plane of the 
sky as measured in 13 CO(2-l), where the disk is least contaminated by surrounding features 
exhibiting non-circular rotation (see §4.2.21 and below). This is significantly different from 
the inclination of the large-scale galactic disk of ~78°. Observations in atomic hydrogen 
(HI) gas also indicate a small change in the inclination of the galactic disk with radius (Ott 
et al. 2001). 



As mentioned in JTJ ICunningham fc Whiteoaki ( 120051 ) measured a size for the central 
molecular concentration in HCN(l-O) of 8'.'5 x4'.'2 (±0'.'5) and major axis at PA = 64°. This is 
only about h alf the size that we measured in 12 CO(2-l) of 16'.'4 x 10'.'8 (±0'/l). As in previous 
observations. ICunningham &: Whiteoaki (120051 ) attribute the observed HNC(l-O) emission to 
the two cross-sections of an edge-on ring or edge-thickened disk. Instead, our results suggest 
that the HNC(l-O) emission originates from the inner region of the disk that we observe here 
in CO. The critical density of molecular hydrogen gas for collisional excitation of HNC(l-O) 
is ~10 4-5 cm~ 3 , which is about an order of magnitude higher than that for 12 CO(2-l) of 
~10 3-4 cm -3 . The excitation temperatures of these two species, however, are comparable. 
The smaller size of the disk in HNC(l-O) compared with that in 12 CO(2-l ) there fore suggests 
a radial decrease in (average) density with radius. iMauersberger et al.l (119961 ) reached the 
same conclusions from a comparison of the line intensities and sizes for the central source 
in 12 CO(1-0), 12 CO(2-l), and 12 CO(3-2). In addition, the average gas density on the north- 
eastern side of the disk is a factor of a few higher than that on the south-western side 

(mm- 



The peaks in HNC(l-O) emission (Fig. 11 of ICunningham fc Whiteoaki 120051 ) span ap- 
proximately the same range of radii as the Pact starburst ring. We therefore associate this 
dense molecular gas with fueling the circumnuclear starburst. The disk that we observe in 
CO extends beyond the starburst ring, and therefore traces more diffuse gas in the disk 
at densities ~10 3 cm" 3 f 33 .2.41) . With an estimated star-formation rate of ~0.4 M yr _1 



( IMoorwood fe Olival 119941 ). compared with a mass in molecular gas for the disk of ~1- 



2 x 10 8 M , the circumnuclear starburst must therefore be a transient phenomenon (lasting 
no longer than ~10 8 -10 9 yrs) if the disk is not replenished from its surroundings. 

The dynamical mass of the disk within a radius of 5" (95 pc), where the disk exhibits 
circular rotation, is ~4 x 10 8 M . Using the method described in §3.2.31 we infer a molecular 
gas mass within this region of ~3.1 x 10 7 M Q (from C 18 0), which is about 13 times lower. 
The dynamical mass within a radius of C'3 (5. 7 pc) of the central s uper massive black hole 
as inferred from H2O masers is 1.0 x 10 6 M dGreenhill et al.lll997h . The central region of 
the disk within a radius of 5" (95 pc) must therefore be dominated in mass by stars or, less 
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likely, by gas not in molecular form. 



4-2.2. A surrounding bar 



The S-shaped asymmetry in the isovelocity contours of the intensity- weighted 12 CO(2-l) 
mean- velocity map (Fig.HJ) rese mbles that s een ex tending to much larger scales in both molec- 
ular and atomic hydrogen gas. lOtt et all (120011 ) have imaged the entire disk of NGC 4945 
in both 12 CO(2-l) and HI at comparable angular resolutions of ~24" with the SEST and 
Australia Telescope Compact Array (ATCA) respectively. They reported a similar anti- 
symmetric distortion in the is ovelocity contour s around the central concentration, extending 
outwards to 100"-200" in HI. lOtt et all (120011 ) attribute these distortions to a (nearly edge- 
on) bar. If correct, then this bar must extend inwards (almost) to the central disk, and may 
be responsible for channeling gas inwards to form this disk. 



4-2.3. The inner kinematically- decoupled component 

The spatially-unresolved but kinematically-decoupled inner component described in 
§3.2.11 spans approximately the same range of velocities as the surrounding disk. We have 
conducted the following check to make sure that this feature is not an artifact of the finite 
angular resolution of our observation. If within a radius of ~5" the central concentration 
can be described by just a disk in rigid-body rotation, then the PV-diagram along its minor 
axis should show only emission at those velocities within the synthesized beam. For a radial 
velocity gradient as measured from the PV-diagrams along the major axis in 13 CO(2-l) and 
C 18 0(2-l) of ~27 km s _1 arcsec -1 , and with a width for the synthesized beam along the ma- 
jor axis of 4'.'4, the PV-diagram along the minor axis should therefore only show velocities 
of 520-640 km s" 1 . Instead, irrespectiveof the position angle of the cut through the center 
of the disk, the inner feature spans velocities of 420-760 km s _1 . 

The gas density of the inner component as inferred in §3.2.41 is sufficiently high to 
excite emission in the J = 1 — > transitions of HCN, HCO + , and HNC as observed by 
Cunningham & Whiteoak (2005). Indeed, their channel map in HNC (which suffers least 
from line absorption against the central continuum source) shows emission at or near the 
centroid of the disk even at the largest blueshifted and redshifted velocities of 370 km s _1 
and 750 km s -1 . Similarly, their PV-diagram in HNC along the major axis of the disk shows 
that, even at disk center, the emission spans a broad range of velocities comparable with the 
range observed for the inner kinematically-decoupled component. Because Cunningham & 



-18- 



Whiteoak (2005) do not show a PV-diagram along the minor axis of the disk, however, we 
cannot be entirely sure that the inner component that we detected in CO also was detected 
in HNC. 

What is the nature of this inner kinematically-decoupled component? The near-IR 
vibrational line of molecular hydrogen traces the walls of a conically-shaped cavity with 
roughly the same lateral si ze as the Pa« starbu rst ring, attributed to a superbubble blown 
by supernova-driven winds ( Marconi et al.lEoOol ). The inner kinematically-decoupled feature 



may therefore correspond to the cooler component of the molecular outflow. If extending to 
the same height as the outflow seen in the near-IR vibrational line of molecular hydrogen 
(~3'/6), such an outflow would have been barely resolved in our observation (with an angular 
resolution of 3'.' 3 perpendicular to the major axis of the disk). As mentioned in §3.2.4} 
we infer a density for this inner kinematically-decoupled component that is 1-2 orders of 
magnitude higher than that of the surrounding disk. This is contrary to expectations if the 
inner component corresponds to molecular gas entrained from the disk or surrounding gas. 

A more attractive explanation is that this inner kinematically-decoupled component 
comprises a separate rotating disk. This may naturally explain the relatively high density of 
~10 5 cm -5 ( §3.2.41) inferred for the inner component. The range in velocities exhibited by the 
inner component is similar to that exh ibited by the water masers around the AGN spanning 



440-800 km s _1 ( IGreenhill et al.l 119971 ). suggesting a connection between the two. If it has 
a radius of 1"4 ( §3.2. ip . the dynamical mass of such an inner disk would be ~2 x 10 8 M Q ; 
if the actual radius is smaller, the corresponding dynamical mass would also be smaller. By 
comparison, the dynamical mass inferred from the H2O mega-masers is 1.0 x 10 6 M , which 
is about two orders of magnitude smaller still. 

The column density of intervening neutral gas derived from X-ray absorption towards 



the AGN is ~10 24,7 cm 2 ( llwasawa et all Il993l ). This is about two orders of magnitude 



higher than the column density of molecular hydrogen gas inferred for the central disk of 
~10 22 8 cm -2 ( §3.2.31) . If the bulk of the abso rption originate s from the large-scale galactic 



disk which in HI has a radius of ~11.4 kpc (jOtt et al.ll200ll ). then the average density in 
this disk is required to be ~140 cm -3 . This is about two orders of magnitude higher than 
the typical density in the interstellar medium of ~1 cm -3 , and so the bulk of the X-ray 
absorption is unlikely to originate from the galactic disk. Instead, this X-ray absorption 
may originate from the inner kinematically-decoupled component. As mentioned in §3.2.4[ 
the denser part of inner component has densities of at least ~5 x 10 4 -1 x 10 5 cm -3 . To 
produce the required column density in X-ray absorption, this part of the inner component 
would then have to span a radius of roughly ~20 pc, and hence have a mass of ~10 M . If 
at this radius the inner component has a rotational velocity of 85 km s -1 , corresponding to 
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half the full range of velocities measured for this component, the enclosed dynamical mass 
within this radius would then be ~10 7 M , comparable to its estimated mass in molecular 
gas. 

The required properties of the denser part of the inner component necessary to produce 
the observed X-ray absorption are comparable with the size and rotational velocity (and 
hence also enclos ed dynamical mas s) of a highly inclined rotating disk in radio recombination 



lines imaged by iRoy et all (120051 ) . This disk has its major axis aligned with the large- 
scale galactic disk, and an ionized gas mass of 10 5 -10 6 M . If the denser part of the inner 
component corresponds to the same disk imaged in radio recombination lines, then this disk 
must be composed primarily of molecular gas. This dense inner component may therefore be 
the hypothesized circumnuclear molecular torus invoked by AGN unification models. In the 
case of NGC 4945, which harbours a Seyfert 2 nucleus, the circumnuclear molecular torus is 
required by AGN unification models to be viewed at a large inclination to its rotation axis, 
as appears to be the case for the disk imaged in radio recombination lines. 



5. Summary and Conclusions 

Previous single-dish maps in the J — 1 — ► 0, J — 2 — > 1, and J = 3 —>■ 2 transitions of 
12 CO infer the presence of a spatially-unresolved molecular ring at the center of NGC 4945. 
The inferred radius of this ring ranges from ~5" at the highest transition to ~15" at the 
lowest transition. We have made an interferometric map that properly resolves for the first 
time the molecular gas as traced in CO at the center of NGC 4945. Our images in the 
J = 2 — 1 transition of 12 CO, 13 CO and C 18 and continuum at 1.3 mm show: 

• A disk in all three lines that exhibits rigid-body rotation within a radius of ~5" 
(~100 pc). Beyond this radius the rotation curve flattens, and the isovelocity con- 
tours exhibit an S-shaped distortion that can be traced in 12 CO(2-l) to a radius of 8'.'2 
(156 pc). There is no central hole in the disk with a radius larger than 2'.' 2 (42 pc). 

• The entire disk has a mass in molecular gas of ~(l-2) xlO 8 M . The dynamical mass 
of the disk within a radius of ~5" (~100 pc), where it exhibits circular rotation, is 
~4 x 10 8 M Q . The mass of molecular gas within this radius is ~3 x 10 7 M , and 



the mass of the centra l supermassive black hole may not be larger than ~1 x 10 M 



(IGreenhill et al.lll997l ). 



Based on the measured line ratios, the bulk of this disk likely has densities ~10 3 cm -3 
and temperatures < 100 K, in agreement with values inferred from single-dish obser- 
vations in multiple molecular species. 
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• The disk extends beyond the circumnuclear starburst ring seen in Paa. The major 
axis of the disk is aligned with that of the starburst ring and also large-scale galactic 
disk, all of which have PA « 45°. The disk is inclined by 62° ± 2° from the plane of the 
sky, significantly smaller than the inclination of the larger-scale galactic disk of ~78°. 

• A newly discovered spatially-unresolved and kinematically-decoupled component at the 
center of the disk in all three lines. This component spans velocities of 420-760 km s -1 , 
which is approximately the same as the range of rotational velocities exhibited by the 
surrounding disk, a smaller disk imaged in radio recombination lines, and an even 
smaller disk in water megamasers. With an upper limit in projected radius as small 
as 1'.'4 (27 pc), this component would not have been recognized in previous single-dish 
maps. 

• The bulk of the molecular gas in the inner kinematically-decoupled component resides 
in two very different densities. The denser portion has a density that is at least one to 
two orders of magnitude higher than that in the bulk of the surrounding disk. The more 
diffuse portion has a density that is comparable to that in the bulk of the surrounding 
disk. 

• A nuclear continuum source dominated by dust with a radius of 4'/ 9 (93 pc), about 
half that measured in 12 CO(2-l). This source extends beyond the bright inner rim of 
the circumnuclear starburst ring to, and along some directions beyond, its detectable 
outer radius. 



By comparison, an interferometric map in HNC(l-O), tracing molecular gas at approximately 
an order of magnitude higher density than the abovementioned CO(2-l) lines, reveals an 
inclined rotating-disk-like feature with a radius of 4" (76 p c), interpreted as a ring with an 
inner radius of 3" (57 pc) (jCunningham &: Whiteoak" 20051 ) . 

We interpret our results in the following manner: 



The primary heating agent of the dust is star formation rather than AGN activity. 
Marconi et al.l (120001 ) infer that UV and soft-X-ray photons cannot penetrate beyond 
a radius of ~l'/5 (~30 pc) from the AGN, whereas the continuum source has a radius 
about four times larger that is comparable with the outermost extension of the starburst 
ring. 

The S-shaped distortion in the isovelocity contours is likely caused by a bar seen nearly 
edge-on. A similar distortion but on larger spatial scales has been se en also in single - 
dish 12 CO(2-l) maps and interferometric HI maps of the entire galaxy (jOtt et al.ll200ll ). 
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In HI, this distortion can be traced out to a radius of ~200" (380 pc), and is attributed 
to a nearly edge-on bar. The S-shaped distortion seen in our CO(2-l) map suggests 
that this bar extends inwards to within ~5" of the center of the galaxy, and may have 
been responsible for channeling gas inwards to create or replenish the central molecular 
disk. 

• The molecular gas mapped in HNC(l-O) traces the inner denser region of the same 
disk that we map in CO. The radial extent of the HNC(l-O) gas is approximately 
the same as that of the circumnuclear starburst ring, and so this relatively dense gas 
is likely responsible for fueling the circumnuclear starburst. In the absence of any 
replenishment, at the present star-formation rate all the available molecular gas in the 
disk would be consumed in 10 8 -10 9 yrs. 

• The denser part of the inner kinematically-decoupled component is a promising candi- 
date for the circumnuclear molecular torus invoked by AGN unification models. If it 
is responsible for producing the bulk of the X-ray absorption seen towards the AGN, 
then this component would have a size comparable with a highly inclined disk imaged 
in radio recombination lines that has a radius of ~20 pc. 

Given the measurements now at hand, it would be instructive to construct a dynam- 
ical model of NGC 4945 to determine whether the inferred large-scale bar is responsi- 
ble for channeling gas inwards to replenish the central molecular disk. As for the inner 
kinematically-decoupled component, we have conducted follow-up observations at higher 
angular-resolutions to clarify the nature of this component, and will report the results in a 
future paper. 

We wish to thank all the SMA personnel in Hawaii, Cambridge, and Taipei for their 
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and S. Matsushita acknowledge support from the National Science Council of Taiwan for 
conducting this work. The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics 
and is funded by the Smithsonian Institution and the Academia Sinica. 
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Telescope & Reference 


12 CO(2-l) J T mb dv 


13 CO(2-l) fT mb dv 


C 18 0(2-l) fT mb dv 




(K km s- 1 ) 


(K km s- 1 ) 


(K km s- 1 ) 


SMA, our results 


925±13 


82±3 


31±3 


SEST, Henkel et al. (1994) 


1050±6 


123.21±3.8 


32±2.0 


SEST, Current et al. (2001) 


740±40 


86±5 


29±2 


SEST, Wang et al. (2004) 


921±0.6 


81.2±0.7 


24.5±0.5 



Table 1: Integrated line intensities measured with the SMA (smoothed to SEST resolution) 
and SEST. 
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Total Flux 


/ T mb dv 


Column Density 


Hydrogen Gas 




(Jy) 


(K kms- 1 ) 


(10 22 cm" 2 ) 


Mass (10 s M ) 


12 co 


1271 ± 20 


713 ± 11 


5.9 ± 0.1 


1.63 ± 0.03 


13 co 


97 ± 4 


111 ± 5 


6.7 ± 0.3 


1.03 ± 0.05 


c 18 o 


29 ± 3 


68 ± 7 


16.5 ± 1.7 


1.2 ± 0.1 



Table 2: Emission intensities and derived hydrogen gas mass for the entire disk 
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Total Flux 


/ T mb dv 


Column Density 


Hydrogen Gas 




(Jy) 


(K kms- 1 ) 


(10 23 cm" 2 ) 


Mass (10 7 M ) 


12 co 


251 ± 5 


2478 ± 49 


2.07 ± 0.04 


3.22 ± 0.06 


13 co 


29 ± 1 


323 ± 15 


1.94 ± 0.09 


3.0 ± 0.1 


c 18 o 


11 ± 1 


128 ± 15 


4.8 ± 0.3 


3.1 ± 0.5 



Table 3: Emission intensities and derived hydrogen gas mass within radius < 5"region 
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Fig. 1.- 

Contours of the continuum emission at 1.3 mm measured with the SMA superposed on a 
negative Paa image of the central region of NGC 4945 taken from Marconi et al. (2000). 
Contour levels are plotted at —3a, —2a, 2a, 3a, and then in steps of 3a to 21a, where a is the 
root-mean-square (rms) noise level of 13 mJy/beam. The main elongated feature in the Paa 
image is the starburst ring. The cross denotes the posi tion of the AGN as determined from 
the position of H 2 megamasers (jGreenhill et al.l 119971 ). The synthesized beam is indicated 
at the lower left corner and has a size of 5'/l x 2'/ 8 and a position angle of 8°. 8. 
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Fig. 2.- 

12 CO(2-l) channel maps with contours plotted from 3a to 75a in steps of 12a, where a = 
0.2 Jy/beam. The velocity indicated in each panel is measured relative to the local standard 
of rest, which is 4.6 km s _1 lower than the heliocentric velocity. The velocity separation 
between each panel is ~20 km s _1 . The cross denotes the location of the AGN, and the 
ellipse marks the approximate size of the starburst ring as shown in the bottom right panel. 
The synthesized beam is shown at the lower left corner of the top left panel, and has a size 
of 5"1 x 2'.'8 and a position angle of 8°. 8. 
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Fig. 3.- 

13 CO(2-l) channel maps with contours plotted from 3a to 24a in steps of 3a, where a = 
0.1 Jy/beam. The synthesized beam is shown at the lower left corner of the top left panel, 
and has a size of 5'.'0 x 2'/ 8 and a position angle of 0°.7. The rest are the same as Figure El 
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Fig. 4.— 

Upper row, left to right: Contours of integrated intensity maps in 12 CO(2-l), 13 CO(2-l) 
and C 18 0(2-l) overlaid on the Paa image from Marconi et al. (2000). Contour levels are 
plotted from 3a with steps of 320a, 32a and 16a for 12 CO(2-l), 13 CO(2-l) and C 18 0(2-l), 
respectively. The ellipse at the lower left corner indicates the synthesized beam. The cross 
denotes the position of the AGN. Lower row, left to right: Corresponding intensity-weighted 
mean-velocity maps in 12 CO(2-l), 13 CO(2-l) and C 18 0(2-l). The velocity measured with 
respect to the local standard of rest is indicated by the horizontal bar at the upper edge of 
each panel. The isovelocity contours are plotted from 400 km s _1 to 800 km s _1 in steps 
of 25 km s _1 . 
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Fig. 5.- 

Top row, left to right: Position-velocity diagram in 12 CO(1-0), 13 CO(2-l), and C 18 0(2-l) 
along the major axis of the disk (PA = 44°). Middle row, left to right: Same, but with a 
PA halfway between the major and minor axis of the disk (PA = 89°). Bottom row, left to 
right: Same, but along the minor axis of the disk (PA = 134°). In all panels, contour levels 
for 12 CO are plotted from 7a in steps of 7a, and contour levels for 13 CO and C 18 from 3a 
in steps of 3a. All velocities quoted here are relative to the local standard of rest, which is 
4.6 km s _1 lower than the heliocentric velocity. 
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Fig. 6- 
Line profile 



in 



12 



CO (2-1) measur ed with the SEST (thin line; from 



Henkel. Whiteoak. fc Mauersbergerl 11994 ) and that measured here with the SMA (thick 



line). The latter corresponds to the spatially- integrated line profile corrected for the primary 
beam of the SMA and convolved to the primary beam of SEST, and using a conversion 
factor of 0.54 to convert from SMA main b eam temperature to SEST antenna temperature 
(IHenkel. Whiteoak. &: Mauersbergerl I1994J ). The root-mean-square (rms) noise level of 
the profiles are 49.6 mK for SEST and 8 mK for the SMA. The integrated line intensity 
measured with the SMA is ~90% (not taking into account any errors in the absolute flux 
calibration) of that measured with the SEST, indicating that we have recovered the bulk if 
not all of the emission present in the same region. 
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Fig. 7.- 

From left to right: Intensity ratio in 12 CO/ 13 CO, 12 CO/C 18 and 13 CO/C 18 along the 
major axis of the disk, measured at intervals of 2" (half the synthesized beam along this 
direction). The error bars indicate ±lcr uncertainty, and the bars with upward pointing 
arrows the 2a lower limits. The negative offset corresponds to the north-east part of the 
disk and the positive offset corresponds to the south-west part. 
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Fig. 8.— 

From left to right: Intensity ratio in 12 CO/ 13 CO, 12 CO/C 18 and 13 CO/C 18 as a function 
of velocity as measured from the PV-diagram along the minor axis of the disk (Fig. [5j bottom 
row). The error bars and the upper limits are the same as in Fig. [7J 
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Fig. 9- 

Results from our LVG calculations for the density and kinetic temperature of the molecular 
gas in the disk. The thin lines are solutions derived from the line ratio 12 CO/ 13 CO for the 
values 5 and 15, encompassing the range of values spanned by this ratio. The dashed lines 
are solutions derived from the line ratio 12 CO/C 18 for the values 13 and 30, encompassing 
the range of values spanned by this ratio. The thick solid lines are solutions derived from the 
line ratio 13 CO/C 18 for the values 2.2 and 3.0, encompassing the range of values spanned 
by this ratio. Below a temperature of ~100 K, the line ratio is only weakly dependent 
on temperature, with smaller line ratios indicating higher densities. Above ~100 K, for a 
given line ratio the density increases with temperature, with smaller line ratios continuing 
to indicate higher densities. 
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Fig. 10.— 

Results from our LVG calculations for the density and kinetic temperature of the molecu- 
lar gas in the inner kinematically-decoupled component. The thin solid lines are solutions 
derived from the line ratio 12 CO/ 13 CO for the values 12 and 17, encompassing the range 
of values spanned by this ratio. The dashed lines are solutions derived from the line ratio 
12 CO/C 18 for the values 17 and 26, encompassing the range of values spanned by this 
ratio. The thick solid lines are solutions derived from the line ratio 13 CO/C 18 for the val- 
ues 1.2 and 1.5, encompassing the range of values spanned by this ratio. The solutions for 
13 CO/C 18 are quite different from those for 12 CO/ 13 CO and 12 CO/C 18 0, indicating that 
both 13 CO and C 18 trace a different denser region than 12 CO. 



